DO-TH-12/26 

Pion production by neutrinos in the delta resonance region and 
possible application to CP searches 

E. A. Paschos* and Dario Schalla^ 
Department of Physics, TU Dortmund, D-44221 Dortmund, Germany 

(Dated: September 19, 2012) 

Abstract 

We present the results of extensive calculations for charged and neutral current reactions of 
neutrinos and antineutrinos in the small Q 2 ^ 0.20 GeV 2 region. The results include single ^7 
and double d Qi dE differential cross sections at energies relevant for oscillation experiments. We 
include nuclear corrections in the Adler-Nussinov-Paschos model and point out that on isoscalar 
43 " nuclear targets there are charge symmetry relations that hold in extended kinematic regions. We 

q; discuss how the results can be used in long baseline experiments in order to study oscillation 

r*j . parameters and search for CP asymmetries. 
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I. INTRODUCTION 



Neutrino interactions have been studied in the resonance region in terms of form factors 
and have been compared with available data [1-8] with considerable success. After the early 
experiments with larger errors, there is a new generation of experiments that provide more 
accurate data. In addition, the oscillation phenomena and the finite mixing angle #13 [9, 10] 
require more accurate theoretical predictions in order to decipher properties of oscillations 
including CP asymmetries. All these require on the theoretical side estimates of amplitudes 
that are reliable. 

One property of production cross sections is the fact that for m^v 3> Q 2 and Q 2 ~ 0(m 2 ) 
the axial contribution is given by the partially conserved axialvector current (PCAC). We 
worked out this framework for resonance production [11] and for coherent scattering on 
nuclei it has been reported by two groups [12-14]. In this article we wish to present our 
results for many reactions, calculate explicitly the energy spectra of the produced pions and 
point out special properties. For example, in the production of the delta resonance in the 
energy region E u = 1.0 to 2.0 GeV the vector squared and the interference contributions are 
almost equal [1 1] . They add up for neutrinos and cancel each other for antineutrinos leaving 
the axial current squared as the dominant contribution in the latter reaction. This will be 
tested in the experiments. It will also be useful for measuring the flux by using this and 
additional channels in the nearby detector and to predict the yield in the far away detector 
in order to establish deviations from normal oscillations, like the presence of a CP violating 
phase. The above property for antineutrino reactions is especially useful because the flux 
of the antineutrino beam will be smaller. For antineutrino interactions we find the cross 
section on a proton target to be smaller than on a neutron. Early qualitative evidence for 
this ordering is already available for integrated cross sections [15, 16]. 

In this article we adopt the neutrino reaction 

u lx +p^fi~ + A ++ (1) 

as benchmark for comparing other reactions of charged and neutral currents. For W < 
1.6 GeV the delta resonance dominates and it has been shown that there is only a small 
background of 10% for other reactions. We shall assume that the I = 3 /2 amplitude dom- 
inates and thus obtain the antineutrino and neutral current reactions. Later on, we may 
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revisit the topic in order to study the changes brought about when we introduce a small 
nonresonant background. 

In section III we calculate the energy and angular spectrum of the pions for the small 
region of Q 2 . For this we introduce the irp — > A — > up cross section in the rest frame of the 
resonance and then transform it to the laboratory frame of the neutrino interaction. Details 
of the calculation are given in section III and in the appendix. 

One worrisome aspect concerns the modifications brought about by the subsequent rescat- 
tering of pions in nuclei. Here we make the observation that charge symmetry predicts similar 
corrections for n + and ir~ interactions on isoscalar nuclei. We shall use this property. For 
nuclear corrections we adopt a model for rescatterings which is both simple and transpar- 
ent [17]. In fact, since some of the neutrino-nucleon reactions are predicted, we can use them 
to test the accuracy of nuclear corrections. 

For isoscalar targets many structure functions are related by charge symmetry at both 
steps of the reaction: the initial neutrino-nucleon scattering and the subsequent pion nucleus 
interactions. Consequently, several relations follow to be tested experimentally. Finally, the 
results for Q 2 < 0.20 GeV 2 are useful for matching them to the form factors at higher values 
of Q 2 . 

The contents of the article are arranged as follows. In section II we define the general 
method and present charged and neutral current cross sections. We point out several regu- 
larities that are inherent in the cross sections. Results for pion spectra in energy E n from 
the axial current are presented in section III. For these calculations we introduce the energy 
and angular dependence of the reaction it + p — > resonance — > tt + p in the rest frame of the 
resonance and transform it to the laboratory frame. The details of the Lorentz transforma- 
tion are given in the appendix. In section IV we compute nuclear corrections for isoscalar 
targets and pay special attention to find charge symmetric reactions. The results will be 
useful for long baseline experiments in which we point out that in our kinematic region 
products of cross sections (g) nuclear corrections will be very similar for the regenerated v e 
and V e beams. Since the mass m e of the electron is very small we repeated the calculations 
with m e and the results are shown in figure 7. The article closes with a general discussion. 

The method can be extended to higher resonances in order to fill the transition region 
between resonances and deep inelastic scattering. In fact, another group [18] uses a dynamic 
hadron model for this purpose. 
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TABLE I: Input quantities and isospin factors for various reactions. 
II. CHARGED AND NEUTRAL CURRENT REACTIONS 

We consider the reactions listed in table I. We take reaction (1) as a standard where for 
X ++ = A ++ there is only the A 3 ^ 2 amplitude. The amplitude for 

Up + n -»■ fjT + A + (2) 

is given by the Clebsch-Gordan coefficient (CGC) shown in the fourth column of table I. The 
CGC is the same for vector and axial-vector form factors and reduces the cross section by a 
factor of y3. When the subsequent decays to specific pion-nucleus channels are considered, 
there are additional CGCs. The results for vector, axial-vector and interference terms were 
computed in an earlier article [11] and the results are shown again in figure 1. Considering 
next the reaction 

77 M + n -> /i + + (3) 

we use the same structure functions as for (1), since they are related by charge symmetry, 
but the sign of the interference term in the cross section changes. We note in figure 1 
that the vector and interference terms are approximately equal and cancel each other for 
antineutrinos, leaving the axial contribution alone. 

In figure 2 we show all charged current reactions on proton and neutron targets. The 
cross section for the reaction in equation (3) is shown in figure 2(d) with its dominant 
contribution coming from the axial current. Figures 2(c) and 2(b) are obtained from 2(a) 
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FIG. 1: Various contributions to the differential cross section at E u = 1 GeV (from reference [11]). 

and 2(d), respectively, by applying CGCs of table I. In all charged current reactions we kept 
the muon mass which is responsible for the turning over of the cross section as Q 2 —> 0. All 
these reactions and especially 2(d) can be used to confirm and also decipher features of the 
reactions. 

For neutral current reactions there are more changes. The effective interaction is 

H eS = (1 - 75) v [xVl + yAl + 7 V°] (4) 

with V 3 and A 3 ^ the isovector and V° the isoscalar hadronic currents. The parameters in the 
hadronic current are given in terms of the weak angle 8\y 

2 

x — 1 — 2 sin 2 9 W , y = —1 and 7 = — - sin 2 6 W (5) 

o 

with sin 2 9w ~ 0.25. The value of y = —1 gives a constructive W3 interference term (because 
of the structure of the lepton current z/7 M (1 — 75) v), making the neutrino reaction larger 
than the antineutrino. Beyond these parameters there is an overall normalization factor in 
the amplitudes. In the charged current interaction appears the current 

Al + l Al = V2(^^)=V2At (6) 

and for the neutral current A^. The CGCs are valid for the triplet (A^ , A^, A~) . Since we 
have taken the amplitude for reaction (1) as the standard amplitude, we divided the neutral 
current CGCs in the fourth column of table I by \/2. 
An additional property of neutral current reactions is 

a(up -> z/A + ) = a{un -> uA°) (7) 
aiVp 77A + ) = aiVn 77A°). (8) 
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FIG. 2: Charged current differential cross sections for E v = lGeV. 

which follows from charge symmetry. The calculated differential cross sections for neutral 
currents are shown in figure 3. The zero mass of the neutrino assures nonzero values for 
the cross sections at Q 2 = which is the exact point from PCAC where neutrino and 
antineutrino cross sections are equal. 

III. ENERGY SPECTRUM AND ANGULAR DEPENDENCE OF PIONS 

Several other variables are measured or have been reported [19-23] in recent experiments. 
Among them are the energy and angle of the pion in the laboratory frame. In our approach 
the lepton part of the interaction factorizes from the pion-target interaction. For this rea- 
son the pion-proton interaction can be computed in the rest-frame of the resonance and 
then transformed to the laboratory frame. For the delta resonance the angular and energy 
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FIG. 3: Differential cross sections for neutral current and E v = 1 GeV. 



dependence is known [24]: 



da 
dVt cn 



a 



(W) 



1 + 3 cos 2 e c 



(9) 



with u v + p (W) given by experimental data [25]. As mentioned already, we are interested in 
the production in the rest frame of the proton, which is obtained by a Lorentz transformation. 
We present the kinematics and the transformation in the appendix. 

Folding this formula with the remaining weak vertex, we derive a triple differential cross 
section 



da^ {E„,Q 2 ,W,E l f>) _G 2 F \V ud \ 



m N 



dQ 2 dWdEf° 



W \EIQ 



L 00 j 2 \ da n + p (W) 



d £lab 



(10) 



The factor before the n + p scattering depends on variables E u ,Q 2 ,u which is typical for 
neutrino interactions. The function Loo defines the square of the weak vertex contributing 
to the helicity zero polarization and includes the mass of the muon. It has been given in the 
articles [12, 13]. We emphasize that formula (10) gives the contribution of the axial current 
alone, which for the antineutrino reaction it is a good approximation, especially in the low 
Q 2 region. For the neutrino induced reaction the vector and interference terms must be 
included as was done in [26], which increases the cross section. 



One obtains 



AE Uh 



from equations (25) to (27) in the appendix by introducing the 



Jacobian which transforms the cross section from the rest frame of the resonance to the 
laboratory frame. 
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FIG. 4: Differential cross section after TU-integration for various values of momentum transfers Q 2 

A (A} 

After integration over W we present the cross section d gf d£lab in figure 4 as a function of 
^lab f or f our va i ues Q f q2 rjj^ dependence on q2 a |. j QW ener gi eS; shown in figure 5, is not 

monotonic but increases up to Q 2 = 0.050 GeV 2 and then decreases. 

Alternatively, we plot d Q2^\^b as a function of Q 2 for various values of E l ^°. In figure 5 
we note that the maximal curve is for E l ^ h = 300 GeV and the peak at Q 2 at smaller values 
than 0.05 GeV 2 . Integrating over E^ h we obtain the curve in figure 6 with the maximum 
region being at Q 2 ~ 0.02 GeV 2 , which a posteriori justifies our small Q 2 approximation. 
The single differential cross section when compared with a recent publication [11] is a 
little smaller. The reason is that the resonant term a n + p (W) of equation (9) is computed in 
this article at W, but the data in reference [11] was computed at cr 7r + p (z/). We repeated the 
calculation evaluating the data as cr n + p (W) and obtained the same result in figure 6 which is 
our prefered calculation. Three comparisons with experimental data [19, 27, 28] were made 
in our earlier publication [11] and more are expected to be done in the future when new 
results become available [22, 23]. 



IV. NUCLEAR CORRECTIONS 



Up to now we presented results for cross sections on free protons and neutrons. The 
experiments, however, are carried out on medium size nuclei inside of which the produced 
resonances decay. Their short propagation in the medium may influence their width [29, 
30]. The main modification, however, occurs in the decay products since the pions can be 
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FIG. 6: Differential cross section after W- and E£ -integration. Note the cange of scale as E v 
increases. 

absorbed by the medium or rescatter exchanging their charges some of the time. We shall 
include these corrections in the Adler-Nussinov-Paschos (ANP) model whose predictions 
have been tested to some extent qualitatively [15]. 

We describe the process in two steps: the production of the pions and then their prop- 
agation through the nucleus. The second step involves a transport matrix that has been 
calculated in the articles [17, 31]. Both steps satisfy charge symmetry. We discuss first 
the results of the model and then we shall remark on the results that follow from the cross 
sections. 

We consider isotopically neutral nuclei where the numbers of protons and neutrons are 
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equal. The production of ir + proceeds through the standard reaction v^p — > fi A ++ — > 
fi~pn + and through 



a{u^n — > fi A + — > fi nir + ) = -a(iy^p — > /i A ++ ). 

9 



(11) 



Thus the total yields on protons plus neutrons in terms of the standard cross section are 



+ 1° / -*++^ 

7T? : 2 -a{u,p /^A ++ ) 
7r,^ : no direct production. 



(12) 
(13) 
(14) 



We have indicated the pions with a subscript % to indicate that they are produced in the 
first step of the interaction. These yields must be folded with the charge exhange matrix 
for C 12 which has been calculated in the original work [17] and also more recently [31]. The 
pions emerging from the nuclei are indicated with the subscript /: 
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(16) 



The overall factor A is for the absorption of the pions and has values from 0.631 for iron to 
0.791 for carbon. 

A similar analysis follows for antineutrinos with the standard cross section being now 
aiv^n — > fi + A~). The structure functions of the two standard cross sections are equal 
through an isospin rotation, however, in the cross section for antineutrinos the vector- 
axialvector interference has an overall minus sign. For incident antineutrinos the yields 
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(17) 



(18) 



We derived the results in the model but several propertiees are more general. The difference 
in the overall cross sections are plotted in figures 2(a) and 2(d). As mentioned already the 
dominant contribution for V^n — > /i + A~ is from the axial current, which in the kinematic 
region of the article was determined by using PCAC. 

We notice that the values in the column matrices in equation (18) are the inverted values 
from those of equation (16). This is a consequence of charge symmetry. It follows now that 



7T 



7T l 



7T 



7T 



(19) 



This relation and other ratios are general. They are valid in other kinematic regions outside 
the PCAC region discussed in this article and they may help the searches for CP asymme- 
tries. Similar results hold for neutral current reactions. We summarize the results for all 
reactions in table II. The results of nuclear corrections given in the tables can be combined 
with equations (16) and (18) in order to give the final yields for each reaction. 



V. ANAYSIS OF THE YIELDS 



The results of this study point to general properties which may be useful to other in- 
vestigations. For the determination of oscillation parameters and the observation of CP 
asymmetries we need the yields of single pions from nuclear targets which have the same 
energy dependence for neutrino and antineutrino reactions. We point out that for isoscalar 
targets the nuclear corrections for tt + and ir~ are equal provided that their energy spec- 
trum at the production vertex are the same. This requires the production by neutrinos and 
antineutrinos to be the same. 
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TABLE II: Yields for neutrinos and antineutrinos including nuclear corrections. 

We pointed out that the structure functions for reactions (1) and (3) are the same by 
charge symmetry. The difference in the cross sections comes from the sign in front of 
W3. We remark now that the contribution of W3 to the cross sections diminishes as E VtV 
increases. This is evident in equation (2.10) of reference [3] and is apparent in figures 6 - 8 of 
reference [26] . The merging together of neutrino and antineutrino cross sections also occurs 
as Q 2 — > 0. We pointed this out for neutral currents in figure 3. It also happens in charged 
current reactions and we demonstrate it for u e p — > e~A ++ and V e n — > e + A~ by repeating 
the calculation and replacing the muon mass by the electron mass. In figure 7 we show 
the two cross sections together and for various energies. For the subsequent interaction the 
nuclear corrections on isoscalar targets will be the same for tt + and 7r~. 

These remarks suggest to use several reactions in the nearby detector in order to determine 
the fluxes and at the same time check the shapes and normalizations of pion production. 
Then measure in the far away detector charge-symmetric reactions. Ratios that follow from 
charge symmetry should provide checks at every step. For instance, in the charged current 
reactions the ir + yield is dominant, but for antineutrinos the ir~ yield dominates. For neutral 
currents the 7T° yield is dominant and the sum ir + + ir~ is, within 5 %, equal to the tt° yield. 

Then one can use the determined fluxes of and in order to compute the regenerated, 
through oscillations, v e and V e fluxes far away. 

For charge symmetric reactions for which the products cross sections ® nuclear corrections 
are the same, the change in the yields at the far away detector should be proportional to 
the change in the fluxes, coming from the oscillations with or without CP violation. 
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FIG. 7: Comparison of muon and electron neutrino or antineutrino cross section for various energies 
E U) u = 1, 2, 5, 10 GeV (bottom to top). The solid line is for electron neutrinos or antineutrinos and 
the dashed one for muon neutrinos and antineutrinos, respectively. 

VI. DISCUSSION 

We carried out an extensive program of calculations in the kinematic region m^v ^ Q 2 
and Q 2 < 0.20 GeV 2 , where PCAC gives the dominant contribution for the axial current. 
We give results for charged and neutral current reactions of neutrinos and antineutrinos. 
One interesting result, in the energy range E u = 1 — 2 GeV, is that the axial contribution 
dominates the V^n — > fi + A~ reaction. In the calculations we kept the mass m M which 
produces the turn over in the differential cross sections as Q 2 —> 0. The long baseline 
oscillation experiments will be searching for u e p — > e _ A ++ type reactions where m e <C m M 
and we repeated the results in figure 7. The PCAC limiting point at Q 2 = is now finite. 

The energy spectra for the pions are also interesting and measurable. In section III 
we outlined a method and carried out the calculations obtaining the energy and angular 
spectra of the pions with the help of transformations summarized in the appendix. Since 
our estimates for the axial current are valid for Q 2 ^ 0.20 GeV 2 we limited the figures to 
that kinematic region. 

For nuclear corrections we adopted the ANP model which is both simple and transparent. 
We presented corrections for the incoherent sum of a proton and neutron target. They can 
be scaled up for other isoscalar nuclei (C 12 , O 16 , ...) by multiplying with the appropriate 
number of particles in the target and using the appropriate transfer matrix. Then we pointed 
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out that some general properties follow from charge symmetry. These relations follow from 
isospin rotations and should hold beyond the kinematic domains introduced by the validity 
of PCAC. 

We hope that the method and results can be adopted in oscillation experiments searching 
for CP asymmetry in the leptonic sector. 
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Appendix: Lorentz transformation from cm to lab frame 

We define as center of mass system the frame where the resonance is at rest. In the 
laboratory the four-momentum of the excited resonance after the collision is 

p'^{u + m N ,\q\). (20) 

We can bring the resonance to rest by the transformation parameter 

P = -J*L_ (21) 
v + rriN 



and the corresponding 



7 = 1 = U + mN . (22) 
7 W 1 ) 



In the rest frame of the resonance with its invariant mass within its width by W, the energy 
of the pion is 

W 2 - m 2 N + ml 

E * = zw • (23) 

With this information we can relate various quantities in the two frames. 
The scattering angles of the pion satisfy 

+ a sin# cm 

7 I cos fl cm + 7-^- 
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The differential cross sections are related by 
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with the Jacobian J derived from the above equation 
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Finally, we obtain 
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